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New semiclassical models of virtual antiparticle pairs are used to compute the 
pair lifetimes, and good agreement with the Heisenberg lifetimes from quantum field the- 
ory (QFT) is found. The modeling method applies to both the electromagnetic and color 
forces. Evaluation of the action integral of potential field fluctuation for each interaction 
potential yields « h/2 for both electromagnetic and color fluctuations, in agreement with 
QFT. Thus each model is a quantized semiclassical representation for such virtual an- 
tiparticle pairs, to good approximation. When the results of the new models and QFT 
are combined, formulae for e and ot s (q ) are derived in terms of only h and c 
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Abstract 

New semiclassical models of virtual antiparticle pairs are used to compute the 
pair Lifetimes, and good agreement with the Heisenberg lifetimes from quantum field the- 
ory (QFT) is found. The modeling method applies to both the electromagnetic and color 
forces. Evaluation of the action integral of potential field fluctuation for each interaction 
potential yields « h/2 for both electromagnetic and color fluctuations, in agreement with 
QFT. Thus each model is a quantized semiclassical representation for such virtual an- 
tiparticle pairs, to good approximation. When the results of the new models and QFT 
are combined, formulae for e and ot s (q) are derived in terms of only h and c. 
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-- 1. INTRODUCTION 

I 

Semiclassical methods have limitations and depend upon classical dynamics, 
which is rigorously valid only for large quantum numbers. Nevertheless, the Bohr model 
and a Bohr-Sommerfeld-type model of meson spectra developed by Brau^ 1 ) exhibit im- 
pressive accuracy even for low quantum numbers. 

It is well known that virtual pairs of elementary particles appear briefly as vac- 
( uum fluctuations and then annihilate each other- 2 ). Many effects of these Virtual An- 
tiparticle Pairs (VAPs) are well understood, such as the phenomenon of screening of 
“bare” electric charge by vacuum polarization of virtual fermions^ 3 ). Can semiclassical 
models of VAPs be constructed? The test of usefulness of such a semiclassical model for 
a VAP should be whether the model successfully predicts the state’s lifetime in agree- 
ment with quantum theory. 


The quantum states of normal material particles obey the Heisenberg Uncer- 
tainty Principle for energy (As At > h/ 2)( 4 ) VAPs are ephemeral quantum states which 
can only exist so long as As At < Tif 2. In general, the maximum lifetime of such a VAP 
(of any massive elementary particle species) is the Heisenberg lifetime 


h/2 

tH = 2^ 


( 1 ) 


where m is the mass of one of the particles (cgs units are used in this paper). 

The goal of this paper is to point out that, if one adheres carefully to the semi- 
classical method, then it can be used to compute VAP lifetimes in good agreement with 


the Heisenberg lifetimes. This shows not only that semiclassical modeling corresponds 
unexpectedly well with QFT, but also offers the advantage that the coupling constants 
appear explicitly in the new resulting formulae, so that we may solve for their values 
(since we know tfj from QFT, too). 

As a tool for further physical insight into VAPs, or a conceptual or teaching 
aid, it would seem natural to apply a semiclassical dynamical model to the VAP (anal- 
ogous to the Bohr hydrogen atom), if it could be appropriately quantized. Semiclassical 
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quantization is a well understood method^ 5 ). The Bohr semiclassical model successfully 
yields so many physical properties of hydrogen and has such undisputed pedagogical 
value that it is still widely used in quantum mechanics texts. In addition, the semiclas- 
sical Bohr- Sommerf eld quantization method continues to produce new physical insights 
into quantum-mechanical problems, such as Brau’s recent success at modeling meson 
spectroscopy using a Bohr-Sommerfeld treatment of constituent quarks^ 1 ). Comparison 
between the prediction of such a dynamical model for VAPs and the Heisenberg lifetime 
(Eq 1) might be instructive in relation to the differences between semiclassical and fully 
quantum-mechanical physics of VAPs. This expectation is encouraged by Brau’s deriva- 
tion of useful formulae for the dependence of some physical quantities on quantum num- 
bers, which emerge from the semiclassical method but have not been obtainable from full 
quantum calculations. Curiously, there seems to be no trace in the literature of such a 
dynamical model of VAPs. 

In this paper it is shown that semiclassical relativistic two-body theory can 
enable us to compute good approximations of the Heisenberg VAP lifetimes tn for the 
massive leptons, the quarks, and electroweak bosons. The agreement of the results with 
QFT is very good, and suggests that the models provide useful insights into VAPs. 

2. VIRTUAL MASSIVE LEPTON ANTIPARTICLE PAIRS 
2.1 Dynamics Corresponding with the Feynman Diagram 

Let us first consider a virtual electron-positron pair (VEPP) unperturbed by 
any external forces, using classical electrodynamics. In second-order QED, processes de- 
scribed by the Feynman diagram in Fig. 1(a) occur An electron with four-momentum 
p, a photon with four-momentum k, and a positron of four-momentum —p — k sponta- 
neously appear from the vacuum at one spacetime point xj = (<x, Xi) and propagate to 
another spacetime point X 2 , where they annihilate; energy and momentum are conserved 
(because the particles are “off mass shell’) 
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We shall take the Feynman diagram in Fig. 1 (a) as the basis for a classical dy- 
namical model, and then ensure that an appropriate quantization condition is met (the 
usual semiclassical approach). 

Of course it is not possible to model Fig 1 (a) perfectly in the classical frame- 
work. In the diagram, the particles and the timelike virtual photon are free, with the 
electromagnetic interaction only acting at the vertices at the top and bottom of the 
diagram. This is the usual QED perturbation- theory approach. In contrast, classical 
physics assumes that the particles move m a potential continuously over time; QED per- 
turbation theory treats the interaction without reference to time dependence^ 8 ). So we 
cannot expect the classical dynamical model to correspond with every feature of the di- 
agram. Rather, the objectives here are to answer the following questions: (1) can such 
a VEPP appear, consistent with the classical principles? ( 2 ) what happens to the par- 
ticles thereafter to determine their dynamics? and (3) is the system subject to reason- 
able quantization conditions as in the semiclassical method? The general viewpoint of 
this approach is that QED has led the way in showing us the occurrence of processes 
like Fig 1 (a) in nature, now that we know about these processes, we can investigate 
what non-perturbative theory might tell us about them - in this case, semiclassical non- 
perturbative theory. 

Concerning question ( 1 ), antiparticle creation and annihilation have been rec- 
onciled with classical relativistic dynamics for a long time. The role of antiparticles in 
classical relativistic dynamics is sometimes misunderstood, because of the historical in- 
terpretation of negative-energy states in relativistic quantum mechanics as antiparticles. 
However, the natural appearance of antiparticles as particles in a state of motion such 
that dr = —dt / 7 (where r is proper time and 7 is the Lorentz factor) was realized by 
Stiickelberg^ 9 ); Costella et al/ 10 ) give a clear exposition of the positive rest mass energy 
of antiparticles in classical xelativistic theory (see also Trump and Schieve^ 11 )) So parti- 
cle/antiparticle propeities are well defined and consistent with the classical principles. In 
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addition, classical-pair creation has been studied, for example by Carati^ 12 ^ in solutions 
of the well-known Abraham-Lorentz-Dirac equation. It should be clear that antiparticle 
pair creation and annihilation are accepted as consistent with classical relativistic dy- 
namics in the way that the following model is developed. 

According to the Feynman diagram, the electron and positron are each created 
at the origin of the center- of-mass (CM) coordinate system. At the exact instant of the 
| particles’ creation, both are at r = 0. In the CM frame the leptons then travel apart 
along a radial line with some initial velocity vq (Fig. lb), with a separation R{t ) = 2r. 
This motion is consistent with the usual theory, since the particles are off mass shell so 
that momentum is conserved. Because the electromagnetic potential is singular at R = 
0, it appears at first as if the pair would be inseparable, due to infinite attraction and 
couldn’t escape from the origin. Let us assume that each particle appears at some finite 
distance from the origin so that their initial separation is R in it, and for the moment let 
us ignore this obstacle. We shall have to let — > 0 to faithfully model the Feynman 

diagram, and we shall find later that R init may be taken to zero in the computations. 

Sometimes the appearance of a VAP is regarded as a violation of the conser- 
vation of energy on a timescale too brief to observe. That is true of the particles’ rest 
mass energy, but since the particles are off mass shell and energy is conserved, the to- 
tal energy of the particle-field system must be zero in QED (it was zero for t < ti). In 
this semiclassical model also, we must consider energy conserved (including rest mass), 
and assume that after the leptons’ emission from the origin, the sum of their rest mass 
energy and kinetic energy is the negative of the electromagnetic binding energy at the 
separation R = 2r of the particles, requiring the two-body system energy equation 

e v = 27m e c 2 + U(R) = 0 (2) 

where 7 = 1/y/l — (v(r, t)/c) 2 is the Lorentz factor (as v — > c, 7 — » 00; a s v -) 

0, 7 — » 1), and U(R) is the potential energy function of the election and positron’s mu- 
tual electromagnetic attraction. Thus our model has zero net energy in the CM frame 
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of reference^ 13 ), and thus it costs the^vacuum no net energy to create these particles. - 
In fact, this is necessary for consistency with the semiclassical approach, which dictates 
conservation of energy in the dynamical model. 

Another way to write eq. (2) makes the difference between VAPs and “real” 
pairs more explicit. As expressed by Greiner^ 4 ); if a virtual pair is separated during the 
interval tjj and if it has gained more energy than 2 me 2 , then the particles become real. 
This is easily exhibited since the usual kinetic energy of an electron is T e = (7 — l)m e c 2 , 
so eq. (2) may be rewritten 

T + U — 2(7 — l)m e c 2 + U = —2 m e c 2 . 

The definition of energy baseline in Eq (2) captures in a classical way the energy debt 
that would have to be paid in order for these leptons to be escalated into free parti- 
cles. The virtual quantum system evolves as if it were in a potential energy well of depth 
2 m e c 2 lower than a “real” pair of particles would be. 

In this model, we shall ignore radiation reaction, since in QED these particles 
are not capable of emitting radiation (and if the particles radiated, then we would be 
modeling some other, more complex process than Fig. 1(a)) We impose conservation of 
angular momentum, so that the spins of the particles are such that they sum to zero net 
angular momentum (the virtual photon is only a carrier of energy, but not of momentum 
or net angular momentum). Then the spin angular momentum of the electron must be 
antiparallel to that of the positron, and the magnetic moments /u e - and /i e + must be 
parallel 

What forces determine the dynamics of these particles? With zero orbital an- 
gular momentum in the system, there are only the electrostatic and magnetic spin-spin 
forces in play (gravity can be ignored). In the rest frame of each lepton, the electro- 
static potential is —e 2 /R and the magnetic spm-spin potential is the magnetic field 
gradient potential energy Ug (see Ref 14), which gives the familiar field gradient force 
V(/i-B) = V( — Ub) (conventionally ignored in discussing the diagram in Fig. 1(a)). 
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- - ' What are the relative magnitudes of the two electromagnetic forces? To' an- 

swer this, we need to estimate the maximum separation vector between the electron and 
positron. Conventionally, an upper bound is calculated for the separation of the electron - 
and positron, which cannot move apart faster than c, the speed of light. Then an upper 
bound on their separation^ 15 ) is cf// e , with f// e given by (recall Eq. 1) 

tHe = T— 2 « 3.2 x 10" 22 s . 

2m e r 

Then ctn& — h/(4m e c) « 10“ 11 cm. 

Let us consider first the case in which the electron and positron separate along 
the direction of ^ e - . The absence of torques ensures that the motion is purely one- 
dimensional. At maximum separation, the potential energy of the spin-spin interaction 
is Ub = — = — 2/ig/j R 3 , where 


ge h 

= 2^2 “ 2eC ‘"' ’ 


( 3 ) 


and g is the Lande factor ~ 2 for fermions. Then the ratio of the potentials is 

Ub_ _ 2^ R / ctfje \ 2 

V e R 3 e 2 ~ V R ) ' 

For R = ct// e , the magnetic force is dominant. Although this differs from the usual 
QED model for VEPPs, it should be noted that this creates no conflict whatever be- 
tween the present model and QED, because diagrams like Fig 1(a) do not contribute to 
any observable effect m QED^ 3 ), other diagrams (processes) in which VAPs interact with 
other particles lead to the observable effects. 

Perhaps it should be noted for comparison that a VEPP is a very different 
bound state of the electron and positron from positronium^ 16 ); an electron and positron 
in a positronium state are in a stationary state with spatial extent of order ao ~ 0.53 x 
10 -8 cm (Bohr radius), whereas the leptons in a VEPP are not. The much smaller scale 
of the VEPP and the 1 / R 3 dependence of the spin-spin force are the reasons that the 
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spin-spin force dominates in this VEPP model, although lt-is-usually a small perturba- 
tion. 

Let us assume that the leading term in the potential is the spin-spin force be- 
tween the leptons, and solve for the dynamics using this term alone; we will then check 
to verify whether this assumption is accurate. 

As the particles separate, the electromagnetic attraction between the leptons 
{ will slow their motions until they reach a maximum separation Rmax where the veloc- 
ity and kinetic energy are zero. While the particles move, we must perform a Lorentz 
transformation from the lepton rest frame to the CM frame, so U(R ) = — 2^/(7 R) 3 
(Jackson^ 17 ) explains how this introduces a change from R to 7 R for this one-dimen- 
sional motion). When the particles come to rest at the turning points, it is easy to solve 
the energy equation (2) because 7=1, and obtain 


eW /3 


4.70 x 10~ 12 cm 


where r e = e 2 /(m e c 2 ) « 2.82 x 10 13 cm is the “classical electron radius”, and A c = 
h/(m e c) « 2.43 x 10 -10 cm is the Compton electron wavelength. 

We see that Rmax ~ c ^He/ 2, which is quite consistent with the conventional 
upper bound on the leptons’ separation derived from QED, and mentioned previously. 

It should be noted that any structure in the pointlike electron or positron is 
known to be smaller than 10 -15 cm, as determined from high-energy scattering^ 18 ), so 
the scale of our model is more than three orders of magnitude larger than any of the 
constituent particles’ structure, justifying our treating the leptons as points. (We now 
explicitly assume < Rmax-) 

At the scale Rmax, the strength of the electrostatic potential V e (R) = —e 2 /R 
is small relative to Ub(R )• Their ratio is Ub/V c ~ 33 at R = R ma x (and greater for 
R < Rmax)- The lepton motion in this model is governed by the spin-spin force, and the 
electrostatic potential is negligible in determining the dynamics. 
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After the leptons reach their turning points, we then expect that the electron 
and positron each retrace their motions in reverse and annihilate one another when they 
collide at the origin. The motions of the particles can be computed by integrating the 
energy equation (2), and we can compare the dynamical result (which we shall call t ve ) 
with the quantum mechanical maximum timescale for this process of creation and anni- 
hilation in the presence of no external forces (the Heisenberg lifetime), tjj e - There is a 
' considerable literature of attempts to solve rigorously the classical problem of colliding 
relativistic charged particles^ 19 ) ; no results are found in the literature for the problem 
considered here. 

Let us also address one conceivable objection to the physical applicability of 
this model. It is well known that a free electron cannot be localized in a wave packet 
smaller than A xz ~ A c /27 r, due to the Zitterbewegung ( 2 °) (oscillations in ( x ) with fre- 
quency u>z = 2 mc 2 /h ~ 1.5 x 10 21 sec _1 ) Although Rrnax <£ A xz, the period of Zit- 
terbewegung oscillations is 4.2 x 10 -21 s i// e , so in QED the Zitterbewegung does not 
have time to delocalize the VEPPs before they annihilate; we shall see below whether 
the same is true in our dynamical model. 


2.2 Solving the Equation of Motion 


Let us calculate the dynamical timescale for the VEPP creation-annihilation 
process, using the energy equation (2) as an equation of motion. We can rewrite the en- 
ergy equation as 



To write the above in dimensionless form, let us define Q = RIR max . Taking the square 
root of both sides of the energy equation, we obtain 



1 - 



Invert both sides and expand 7: 
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Solve for v: 



an equation of motion which may be integrated. 

Let the total time for the formation, separation and re-annihilation of the vir- 
tual e + e _ pair be denoted t ve , the dynamical lifetime. Then t ve is twice the time for the 
particles to separate and reach the turning points, since the path is approximately time- 
symmetric. For the purpose of this calculation, we can let R zni t — > 0 without encoun- 
tering a singularity, and the model now corresponds to the Feynman diagram without 
depending upon a free parameter i£* m t . Then 



The integral may be found in ref. 21; substituting its value yields the dynamical VEPP 
lifetime 


\/^ r (|) 

6r d) 


-(1.7247) = 


e 2 h 2 \ 1/3 1.7247 


2.70 x 1CT 22 s . 


Let us compare this dynamical lifetime with the Heisenberg lifetime for such 
an electron-positron pair, t// e = h/(4m e c 2 ) « 32 x 10 -22 s above. The dynamical 
computation is 16% lower. Since t ve « i// e , there is clearly no difficulty for the model 
posed by Zitterbewegung, in agreement with QED. 

We were trying to obtain the maximum lifetime of such a dynamical system, 
and the reader may wonder whether t ve is in fact a maximum If the electron and posi- 
tron separated along some trajectory different from the torque-free axis of their angular 
momenta (which we used in the above calculation), what would the lifetime be ? Let us 
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refer -to Fig. 2. Let us suppose that the electron and posit-ron initially separate along 
a line with a direction at an angle 9 relative to the /u vector. If 9 = 90°, then Ub = 

— fj. B > 0 and the energy equation (2) does not have solutions because the potential 
is repulsive, not attractive. The semiclassical model suggests that VEPPs do not exist 
with motions along these trajectories. The cones such that 9 & 54.7° define the sur- 
face outside which Ub > 0, as shown in Fig. 2. Inside the cones, as 9 decreases n B 
becomes positive (Ub becomes negative) and it becomes possible to satisfy Eq. (2) iff 
/a - B > 2 m e c 2 . 

The dynamics of the VEPP are more complex for 0 < 9 < 54.7° than in the 
9 = 0 case, which we already treated. Let us consider an electron near the vertical axis, 
as pictured in Fig. 2, such that 0 < 9 <g. 54.7°. There is now a torque on the magnetic 
dipole moment of the electron, given by r = ^txB. Classically this causes the magnetic 
moment to precess about B. But by symmetry (and conservation of angular momentum) 
the positron experiences the same torque due to the electron’s magnetic field, and hence 
both magnetic moments precess equally and remain always parallel. Because n and B 
are no longer parallel, |/i-B| is reduced. If the trajectories of the electron and positron 
are still close to the vertical axis, R in the particle frame still is approximately 7 R in the 
CM frame, and we can use Ub = — / i 2 ( 3 cos 2 0 — l)/(7 R) 3 in Eq (2)( 22 ) As 9 increases 
from 0 to larger values, this leads to smaller values of Rmax and t ve , and thus we see 
that the values for these parameters computed first above are indeed maximal as desired. 

We see that the conservation laws uniquely determine the dynamical model 
that yields t ve classically, giving us a new derivation of tfj e to good approximation. It 
will be demonstrated below that the agreement of t ve and t^e is no coincidence, but re- 
sults from valid physical relationships. 

Let the numerical factor involving the ratio of F functions in eq. (5) be desig- 
nated f s & 1 72. Note that for the ratio 

tf-Ie _ I* 1 / 4c \ 

T v 7 " 4j s \^h*J 


« 1.16 
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m e cancels out, so this form of VAP lifetime derivation works for all of the massive lep- 
ton VAPs (electron-positron, muon-anti-muon and tau-anti-tau leptons), irrespective of 
their-masses (as it must to be valid). The muon and tau lepton, of course, have smaller 
values of R m ax and t v because of their larger masses, and t vil ~ t VT ~ t// T to the 
same relative precision. (See Table I.) 

These results show that this semiclassical model gives remarkable agreement 
with the QFT lifetimes of VEPPs, in view of the model’s simplicity (only point magnetic 
dipole particles interacting via a potencial function in one dimension). But the reader 
may wonder whether this is coincidental. And what about the quantization that is nec- 
essary to make it really “semiclassical” rather than entirely classical? Next the quantiza- 
tion will be addressed, and an independent physical derivation of the VEPP quantization 
will be presented. 

2.3 Quantization Condition is Met 

The Bohr model was quantized using conditions on the angular momentum of 
the orbiting electrons in a non-relativistic case^ 23 ). In the present model, however, the 
orbits have vanishing angular momentum and angular momentum conservation has al- 
ready been imposed to determine the orientations of the n e - and n e + . Hence the Bohr 
approach is not exactly applicable in this VEPP model. 

Let us consider the quantization of the electromagnetic field in the VEPP. The 
well-known zero-point field energy of the electromagnetic field is e ern = ^hu. That is, for 
any wave mode of angular frequency u>, there is a vacuum fluctuation energy £ ern . An- 
other way to view this fact is that there is a quantization of any zero-point fluctuation 
mode oj such that its energy is proportional to the action J h . (Quantization of action 
is a fundamental connection from classical mechanics to quantum mechanics, and was 
-actually used to derive the Schrodinger equation from the principle of least action by 
Feynman^ 24 ).) Since the VEPP induces an electromagnetic vacuum fluctuation, let us 
calculate the action associated with the electromagnetic field in the VEPP over its cycle 
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from creation to annihilation 

The expression for the action associated with a potential U acting on a particle 
is given (for the relativistic case) by Lanczos^ 25 ) 



where ds — c dt/j. It follows that the action integral of the VEPP electromagnetic field 
' is (recall the time symmetry relative to the turning point at t ve /2): 



4^2 R r 


c R z 

C - £t mox 


*cVi - e 3/2 c R2 m 


f 1 dC_ 

Jo C _3 C 3 \/l 


c- 3 c 3 vi - c 3/2 


2 m e c fi\ 

R max m eC 2 


f'- 

J o Ji 


r, e o 2.* f ^ ^ 

= 2 m e c fa — 2 m e c t ve = 2 

_ £3/2 C J VB \ 4 C 


2*2 \ 1/3 


This final result 2 is approximately 4.4 x 10 -28 erg s, 16% less than \Ti Thus the model 
is inherently quantized (to within the 16% difference in t ve ), and is the corresponding 
semiclassical model for VEPPs. 

In the following sections it will be seen that when the same semiclassical 
method of modeling is applied to quark- antiquark virtual pairs, the same agreement in 


the quantum mechanical and semiclassical timescales results, as well as the action inte- 


2 This action integral for our model VEPP potential equals the conventional zero- 

point field action, to good approximation. N.B: if t ve exactly equalled t// e , then A would 
exactly equal j h . In Section 6. below an interpretation of this is given. 
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3. VIRTUAL QUARK ANTIPARTICLE PAIRS 

3.1 Dynamics Corresponding with the Feynman Diagram 

The method of Sec. 2 can be applied to quark-antiquark virtual pairs A pro- 
cess analogous to Fig 1(a) exists in Quantum Chromodynamics (QCD), but with a 
gluon replacing the timelike photon^ 26 ). The energy equation is of the same form, 

e q = 2'jm q c 2 + U(R) = 0 , (2') 

except that the dominant force is the “color” force instead of the spin-spin electromag- 
netic attraction. The potential function TJ ( R ) has been inferred from decades of high- 
energy scattering experiments. Mesons like the J /ip particle are bound states of a quark 
and its anti-quark, and high-energy particle collision experiments have established phe- 
nomenological forms of U(R). These experiments have revealed that a useful empirical 
form of the potential is the so-called “Cornell funnel potential” ( 27 ) 

V(R) = ~ + aR, 

where k ss |a s , and the QCD strong coupling strength ~ 0.2 in “natural units” (units 
such that h = c = 1), and a « 0.25 GeV 2 . Converted to cgs units used herein, k ~ 8.6 x 
10 -18 erg cm. The second term, aR, is only significant for R > 10“ 13 cm, so we will not 
need to consider it, since the first term with the Coulomb- like dependence will strongly 
dominate the potential for the range of R mox of quark VAPs. The Coulomb-hke term 
is due to exchange of a vector gluon between the quark and antiquark^ 27,28 ). (There are 
other empirical forms of the quark-antiquark potential - many discussed in Ref 27 - but 
the Cornell potential will be sufficient for present use. Also, it should be mentioned that 
the light up and down quarks are not as accurately treated with this kind of potential 
model as the heavier quarks; nevertheless, it is well known that potential models work 
surprisingly well in describing the dynamics of all but the lightest quark flavors.) 
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The QCD coupling coefficient a s is not constant, but is a “running coupling” 
because of its variation with 4-momentum transfer q of the interacting particles^ 29 ) . This 
can be allowed for in our semiclassical model However, Jt is important to be clear about 
what this two-body model can be expected reasonably to correspond with in QCD The 
QCD coupling strength a s {q) has been calculated to four- loop order ( 3 °), including a 
series of processes (diagrams) with numerous topologies. But the present semiclassi- 
i cal model only consists of the spacetime structure of a single-loop process. Hence it is 
not applicable to processes in which more than one loop are involved. It makes sense to 
compare the two-body model here only with lowest-order QCD (and corrections with the 
same spacetime diagram, such as spin-spin or electrostatic corrections). By the same to- 
ken, in QFT a quark VAP doesn’t include those multiple-loop processes, so QFT logic 
prohibits using multiple-loop processes to compute tfj q for quarks also. 

First, to demonstrate the basic physical model, let us take a s = 0.20, some- 
times called the “canonical” valued 27 ), derived from a fit of the charmonium system. Let 
us calculate the dynamical timescale for a charm quark- antiquark VAP using this single 
value of a s . After that, we will then generalize the model to account for the running of 
a s as a function of q. 

Here is the semiclassical dynamical timescale calculation: including a Lorentz 
factor 7 appropriate for the 1/R color-Coulomb function (to make the Lorentz transfor- 
mation to the CM framed 17 )), the potential in Eq (2') becomes 


£/(*) = 


K 

7 R 


As in Sec. 2, we can solve (2') for Rmax a t the turning point 
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the color-Coulomb part of the potential, and the aR term of the- Cornell potential is in- 
significant in magnitude (as is the electrostatic potential). The dynamics are determined 
by the single vector gluon exchange, not the multi-gluon exchange oTconfinement m the 
linear region, in agreement with QCD^ 28 ). 

For the dynamical timescale, we rewrite the energy equation (2') as 7 2 = , 

with C = R/Rmax > and find 

Rmax/2 dr __ Rmax f 1 d£ 

c^W c Jo V'W 

_ Rrnax _ R TTL CL 2C p f° K (Zl\ 

~ c r(§) " C Jc ~2 m q c* 1 ' 

where we label the £ integral’s value as f c « 2.0. For the charm quark, with m q « 2.7 x 
10 -24 g, we t vq « 1.3 x 10 -25 s. 

For comparison, the Heisenberg lifetime of the charm quark is £// q ~ 1.1 x 10~ 25 
s, about 14% lower. 

Considering the quantization of the color force potential field in a virtual charm 
quark-antiquark pair, we compute the field action over the VAP lifetime: 

W 2 / k \ c dt _ 2k 2 dr 

V C 1 R ) 7 c io 7 2 i?^l - C 

k r 1 d^ _ k c 1 d£ 

~~ c Jo c - 1 C\/w = c Jo yr^c 

K/ 2/c 

= -f c « — » 5.6 x 10~ 28 erg s « 0.53 x % ( 6 ') 

This action integral is a good approximation to the zero-point fluctuation action, - h , 
and means that this version of the model, too, is inherently quantized just as the elec- 
tromagnetic case was in Sec. 2. N.B: if t vq exactly equalled t// q , then A would exactly 
equal ^h. 
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3.2 Generalization to Running Coupling a s (q ) 

The above results illustrate the basic application of the model to a quark VAP 
for which a s was determined from experiments with real (non- virtual) bound states of 
that quark. But to make the model realistic for quarks in general, the running coupling 
a s (g) must be included in the model, not just one universal value of a s We can develop 
the model to account for- the running coupling if we can include a characteristic q for the 
VAP of each quark flavor 

The QCD running coupling is a function of momentum transfer Q 2 = — q 2 . A 
standard way to incorporate asymptotic freedom is to let a s depend on the quark masses 
and use — q 2 = (mi+m 2 ) 2 , with the m t the quark masses^ 28 ). Then we may use Richard- 
son’s expression^ 31 ) 

33-2n / ln(-g 2 /A 2 ) ’ 

where nf is the number of quark flavors with masses much less than m\ + m 2 , and A is 
the QCD scale energy, approximately 0.22 GeV. 

The parameter nj varies as described above because we are using an effective 
field theory obtained by integrating out all quarks heavier than the typical energy under 
consideration^ 32 ) I.e., rif is the integer such that there are nj quarks with masses much 
less than 2 m q ; for this paper nj will be defined as the number of quarks with masses of 
2m q /10 or less. We now can use mi = m 2 = m q in the expression for — q 2 and com- 
pute k from eq. (7) more realistically for a running coupling. The results for this run- 
ning color-Coulomb force model, with a a (q ) given by Eq. (7), appear in Table I (marked 
with the symbol © to indicate “Color Only”). 

We see that the semiclassical lifetimes t vq of the quark VAPs all agree with 
the Heisenberg lifetimes tn q from QFT to within a factor of about 2.8, and much bet- 
ter for the charm and bottom quarks. While this agreement isn’t as good as for the lep- 
ton VAPs in Sec. 2, the semiclassical model is still in rough agreement with QCD, to 
much better than an order of magnitude. The color-Coulomb potential from QCD, when 
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taken as a classical force strength, leads to dynamical timescales of VAPs that 'corre- 
spond surprisingly well with the QFT lifetimes for the VAPs. In addition, the values of 
R^nax are consistent with VAP size scales that are deep within the color-Coulomb part 
of the potential function (not out in the linear region where the theory would be inap- 
propriate). This needn’t be a total surprise in view of Brau’s^ recent success modeling 
mesons semiclassically, but the present work is the first semiclassical treatment of quark 
VAPs. 


S.3 Combined Color-Coulomb and Spin-spin Interaction Models 


The modeling of quark VAP lifetimes can be taken as a qualified success, but 
it is possible to go further with this model and include the effects of spin-spin interac- 
tions between the quark- antiquark pairs As noted by Lichtenberg^ 33 ) , the spm-spin in- 
teraction between quarks is the source of electromagnetic mass splittings among hadron 
isospin multiplets. In the present semiclassical model, because the quark VAPs have 
sufficiently small Rmax values, the spin-spin energy can be comparable with the color- 
Coulomb potential for the heavier quarks in which asymptotic freedom reduces a s (q). 
Let us include the spin-spin potential in the quark VAP energy Eq. (2') and examine 
the effect on t vq 

K 2/ig 

7 R ( , yR) 3 ’ 



where 


_ gQg ft 

^ q 2 m q c 2 


is the magnetic moment of a quark with charge Q q . The energy equation becomes 


e q = 27 m q c 2 


K 


2/i 


9 _ 


= 0 . 


( 2 ") 


7 R (7 R) z 

Again we can derive the turning point separation Rmax by letting 7 — » 1. This is easily 
done by defining £ in terms of the scale length from the previous case which considered 
only the color-Coulomb force. 


C = 


R 

k /(2 m q c 2 ) 





The only difference in the limits of integration from the color- Coulomb-only case is that 
the upper limit becomes Cmoi > 1 instead of just 1. The dynamical timescale then be- 



comes 
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This expression reduces to the previous one (Eq. (5')) when A = 0. 

For each of the six quarks, one computes A(a s ), the value of the dimensionless 
turning point separation of the quarks £ max from the root of Eq. (4"), and then t vq is 
given by the above expression with £ max as the upper limit of integration. The integral 
is simple to perform numerically (a short C program for this is available from the author 
upon request). The results are tabulated m Table I, next to the tfj q values. 

It can be seen that the spin-spin interaction has insignificant effect on t vq for 
the light down, up and strange quarks. Because of asymptotic freedom the QCD cou- 
pling weakens for the heavier quarks, however, and the spin-spin correction comes into 
play, increasing Rmax over what it would be in the case of the color- Coulomb poten- 
tial alone. The top quark VAP, with A « 2.0, is dominated by the spin-spin potential 
rather than color- Coulomb. The net effect of including both potentials is to improve 
the agreement of t vq with tfjq significantly for the top quark, but increase the discrep- 
ancy slightly between the model and Heisenberg lifetimes for the charm and bottom 
VAPs. The model lifetimes t vq for charm, bottom, and top quark VAPs differ from the 
Heisenberg values by 37%, 7%, and 9%, respectively. Including spin-spm attraction has 
insignificant effect on the model lifetimes for the down, up, and strange quarks, which 
differ by factors of 2.8, 2 8, and 2.0 respectively from the Heisenberg values. Potential 
models for these light quarks are well-known to be generally poorer than for the heavy 
quarks ( 27 ). 

4. VIRTUAL WEAK BOSONS 

The third class of elementary particles to which we can apply this model are 
the weak interaction gauge bosons W ± and Z°. These particles are believed to be non- 
composite and pomt-like. It would be interesting if we could use the virtual weak boson- 
antiboson pairs to probe the weak interaction, but our model requires that we use as 
the potential the strongest force to which the particles respond; in this case, the weak 
bosons interact by the electromagnetic spin-spin force as well as the weak force (which 
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is of order ] O 30 times weaker than the electrostatic attraction of the antiparticles in low- 
energy interactions^ 34 )) Consequently the potential that we must use is the spin-spin 
interaction between two massive bosons via their spin magnetic moments: 


Hw = 


2 TTl\yC 


h, 


( 3 ') 


where the factor g « 1 can be assumed for bosons. 

The form of Eq. (3') for the magnetic moment is therefore the same as Eq (3), 
except for the substitution of the W or Z mass for a lepton’s mass, and so the subse- 
quent mathematics is the same as in Sec. 2. (We may wonder whether the Z° has a 
magnetic moment like the electrically neutral neutron, but if so, then the present model 
would seem to be applicable.) We find that the dynamical timescale t v yy is again 16% 
less than tuw- The computation of the action is also 16% less than \ti. Therefore, this 
third type of elementary point-like massive particles also can be described approximately 
by this semiclassical model. (Results are in Table I.) 


5. COMPOSITE PARTICLES AND THE BREAKDOWN OF THE MODEL 

The models above were built around the dynamics of only two point-like bodies 
interacting through a simple potential. Consequently the models should not be expected 
to represent well the dynamics of composite particles like mesons and other hadrons. For 
instance, virtual pairs of anti-mesons comprise four quarks m toto, and this is no longer 
a two-body problem. Nucleons exhibit highly complex internal dynamics due to their 
comprising numerous valence quarks and virtual strange quarks^ 35 ), so our models aren’t 
applicable. If the two-body models that have been constructed above do not describe 
such multi-body VAPs, then that is no objection to the two-body models. 


6. THE COUPLING CONSTANTS 


Having found that the dynamical timescales in these models approximate well 
the Heisenberg lifetimes in these eleven cases - indeed, in all of the applicable cases of 



Virtual Antiparticle Pairs 2.3 David Batchelor July 27, 2001 

known massive non-composite VAPs in nature - let us proceed with some reasonable 
derivations based on the observations t// ~ t v for each case. 

6.1 The Unit of Charge e 


The value of the action integral A in the VEPP case (Sec. 2) implies that there 
is not just a fortuitous coincidence between the tj/ e and t ve timescales. Rather, the 
VEPP model matches the expected action integral for a vacuum fluctuation, and this 
is no doubt because it reflects the physical process in a valid way. This agrees with stan- 
dard QED, in which vacuum fluctuations of the positron-electron field generated by the 
electromagnetic field are the virtual creation and annihilation of VEPPs, and the fluctu- 
ations of the electromagnetic field generated by the VEPPs are the virtual emission and 
absorption of photons, as Schwinger exp lamed ( 36 ). 

Since we know the equation for tfje from QFT, we may equate it with the semi- 
classical result in Eq. (5) to derive approximations for two familiar fundamental physical 
constants, the electric charge quantum e and the fine structure constant a = e 2 /(he), in 
terms of h and c only - a very interesting result. If we approximately equate tn e ~ t ve , 
label the charge of a virtual positron e v to denote the derived value, and solve for it, 
then we find easily that 

he \ 1 ^ 2 

16 / 3 ) ~ ^ ^ X 10 6SU 



This is 30% high, compared with the usual measured value, e « 4.8 x 10“ 10 esu. It is 
closer to the “bare” charge - what e would be if it were not “dressed” by the screening 
of vacuum fluctuations^ 37,38 ) . In fact, we should have expected that this calculation was 
likely to yield a result larger than the usual observable charge e because Rmax -C A c , the 
characteristic scale of the screening cloud of VEPPs which surrounds any free electron 
or positron. Different calculating methods lead to different estimates for the magnitude 
of the bare charge: ref. 37 gives a calculation suggesting that e 2 ore « 1.08 e 2 , but ref. 38 
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presents an argument based on the photon propagator, showing that the QED pertur- 
bation theory fails and e < 1are could be infinite. If the semiclassical model in the present 
paper had led to a value for e v less than the observable 4.8 x 10 -10 esu, then the model 
would appear unphysical; however it is unclear whether the semiclassical approximation 
should be expected to hold to such precision. The value for e v derived from this model is 
consistent with e<, are from QED, given present theoretical uncertainties. 

We can also easily derive a v fn 1/(16/^) ~ 1/83; this is 66% larger than the 
well-known 1/137 037 value for a (just because of the propagation of the 30% difference 
between e v and e). As pointed out in the preceding paragraph, our model tends to re- 
flect the bare value of a, not the usual observed value based on e for distances from a 
free electron that are large relative to A c . 

6.2 The QCD Coupling a s (g) 

The semiclassical model is successful at predicting tfj q within 10% or bet- 
ter for the bottom and top quarks, so it is obviously possible to reverse the derivation, 
start from t vq R 2 tfj q and derive a s (q ) in the one-loop approximation of Eq. (7) for 
\q\ = 2 m q c in the cases of the charm, bottom, and top quarks (that is, |g| = 3.0, 9.4 
and 348 GeV/c). It is straightforward to pick a trial value of a s , calculate A and Cmaxi 
perform the numerical integral to compute t vq accounting for both color and spin-spin 
potential terms, and iterate to a convergence of the predicted t vq with tn q . 

The prediction for the top quark of a: s (348 GeV/c) ~ 0.14 is an interesting one 
from semiclassical theory, in an energy domain not yet reached by particle accelerators; 
it’s about 27% greater than the 0.11 that Richardson’s expression (eq. 7) predicts. For 
the bottom quark, the model implies a s (9.4 GeV/c)?s 0 18, about 5% less than eq. (7). 
For the charm quark, the derived value of a a (3 GeV/c)« 0.15 is approximately 32% less 
than the Richardson expression result, 0.22. 

Although values of a s that are derived at the light quark masses are less ac- 
curate, they are still much better than ordei of magnitude. Since A is small for the up, 



Virtual Antiparticle Pairs 


25 


David Batchelor July 27, 2001 


down and strange quarks, the color-only formulae apply. To derive a value of a s from 
the model, eq. (5') with f c = 2.0 is set equal to t^ q . This gives a s « 2.4x 10 16 m (? c 3 t// 9 « 
0.19 For the strange quark, this is 50% less than Richardson’s value (eq. 7). For the up 
and down quarks, the model implies or 3 « 37% of the Richardson value. 

6.3 The Electroweak Bosons’ Coupling 

A separate weak interaction force constant cannot be derived from the model, 
even though we know t v w ~ t^w- In these cases, though, Rmax 10 -14 cm, the elec- 
troweak unification scale^ 39 ). The electroweak unification in the Standard Model at this 
scale means that it is reasonable for the dynamical model to behave the same as in the 
electromagnetic dynamical model of VEPPs and other massive lepton VAPs. The result 
here is then consistent with the Standard Model, given electroweak unification. 

7. DISCUSSION 

Let us summarize the results. 

(1) The dynamical timescale t v predicted from VAP motions due to the electro- 
magnetic spin-spin force approximately equals the Heisenberg lifetime of the 
VAP, within 16%, for the massive leptons. 

(2) The action integral of the electromagnetic spin-spin potential over the dynami- 
cal timescale approximately equals the zero-point field action, Ti/2, within 16%, 
for the massive leptons 

(3) The dynamical timescale t v predicted from VAP motions due to the quark- 
antiquark color- Coulomb and spin- spin forces approximately equals the Heisen- 
berg lifetime of the VAP tn, within 10% or better for the heaviest quarks (bot- 
tom and top), and within a factor of about 2.8 or better for the lighter quarks 
(down, up, strange and charm). 

(4) The action integral of the quaxk-antiquark color potential over the dynamical 
timescale approximately equals the zero-point field action h/2, within 6%, for 
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the charm quark 

(5) The dynamical timescale t v predicted from VAP motions due to the electro- 
magnetic spin-spin force approximately equals the Heisenberg lifetime of the 
VAP tH, within 16%, for the weak bosons. 

(6) The action integral of the electromagnetic spin-spin potential over the dynami- 
cal timescale approximately equals the zero-point field action h/2, within 16%, 
for the weak bosons. 

Several points should be noted about the above results: 

• There are no free, adjustable parameters in any of the models. 

• (2) does not follow trivially from (1) - i.e., the close agreement between A and 
h/2 has a different mathematical origin than the agreement of t v with tfj. 

• (4) does not follow trivially from (3), similarly. 

• (3) and (4) are independent of (1) and (2) because the color force potential is 
completely different in form from the spin-spin potential. 

• (5) and (6) do not follow trivially from (1) and (2), because bosons and 
fermions have different spins which only lead to similar expressions for their 
magnetic moments incidentally. 

There appears to be no way that these t// and t v values could all be so nearly 
equal by coincidence alone, much less for the field quantization conditions to magically - 
come out satisfied by chance. In the case of the VEPP model, some algebraic depen- 
dence on h exists in the formulae, but not in such a way as to account for the satisfac- 
tion of the quantization condition. The agreement between A and h/2 suggests an im- 
portant relationship between classical and quantum aspects of action. In the case of the 
quark-antiquark model, however, the color potential function is a purely empirical for- 
mula, with no algebraic relationship to h, consequently, the agreement of A and h/2 for 
the charm quark VAP could only be a perverse numerical coincidence if it did not sug- 
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gest that the model has physical significance 

The approximate equality of tu and t v for so many particle types suggests 
that in a fully quantum mechanical theory, the two would indeed be identical. If tjj = 
h/(4mc 2 ) = t v exactly, then the VAP potential field action A = exactly (compare 
the definition of f// and eqs. (6) and (6')). The near-equalities suggest that the mod- 
els represent well the dynamics of wave packets in VAPs, just as the Bohr-Sommerfeld 
1 quantization represents well the dynamics of electronic wave packets in hydrogen atoms 
and constituent quarks in mesons^ 1 ). Of course quantum mechanics and QFT cannot be 
equalled by semiclassical methods, but all of these methods have distinguished roles in 
the understanding of physical problems. 

It is quite intriguing that the quantum fluctuation action \ h emerges naturally 
from these model calculations, because source-free considerations of the equations for 
the electromagnetic field are generally used to obtain the zero-point energy in each wave 
mode, ^hur, in contrast, in the present derivations, the particles are the zero-point field 
sources, as Schwinger showed in QED^ 36 ). 

It is also interesting that we have been able to let R tn %t — > 0 in our compu- 
tations, but it apparently does not make sense for Rmit to be smaller than the Planck 
length (hG/c 3 ) 1 / 2 « 1.6 x 10 -33 cm. At this scale, the smoothness of spacetime 
is thought to break down into a foam-like topology because of quantum gravitational 
fluctuations^ 40 ). This lower limit has no significant effect on the t v values because i? mai 
is much greater than the Planck length in every case. (The gravitational force is not 
amenable to a semiclassical model of the type that has been developed herein, because 
the VAP models have the symmetry property that antiparticles attract but like particles 
repel; gravitation is different in that, like particles attract as well as antiparticles. If one 
writes an equation like the energy equation (2) for the gravitational “potential”, then 
the Rmax is less than the Schwarzchild radius of the VAP masses, indicating that the 
model is not physical ) 
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One possible improvement on this work is to reconstruct the models along the 
more elaborate lines of Oliver’s detailed self-force treatment for the electron^ 41 ). How- 
ever, that would lead to mathematical developments beyond the scope of the present 
paper. Another approach is to consider how the Dirac two-body equation might be 
adapted to this problem. Also, in the cases of quark VAPs, it may be that simplified 
equations for the QCD interactions, which have been studied by Wilczek, could be appli- 
cable as alternatives to the Cornell potential; and improve the poor agreement between 
the model t v of the light quark VAPs and their Heisenberg lifetimes. 

Because the force constants e and k can be expressed in terms of h and c, the 
question arises: why do the electromagnetic and color forces have the strengths that 
they have? Are other strengths conceivable for these forces in some hypothetical alter- 
native universes, as has been suggested by proponents of the “anthropic principle?” ( 42 ) 

It is risky to generalize from only two instances, but these semiclassical VAP 
models for electromagnetic and color force interactions have a clear common govern- 
ing principle: the dynamical VAP timescale must equal the Heisenberg VAP lifetime. 

We may infer that both the electromagnetic and color force constants take on their re- 
spective values because of the i// ~ t v constraints, which imply unique values for e and 
ct 3 (q). This suggests that ft and c are the more fundamental quantities in nature. 

This semiclassical modeling method clearly gives correct answers, to good ap- 
proximation, for the two-body VAPs to which it is applicable. The method therefore 
deserves to be part of the toolkit for understanding quantum theory of virtual particles. 

The author gratefully acknowledges helpful discussions of this work with C. J. 
Crannell, J. M. Hollis, R G. Lyon, R J. Drachman, M. L. Rilee, L. N. Foster, and 
I. T. Durham. Thanks to J. S. Langton for detecting a computational error in an ear- 
lier draft, and for improving the C program used for numerical integration in Sec. 3.3. 
Constructive recommendations from two referees are also greatly appreciated. 
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FIGURE CAPTIONS 

Figure 1. (a) Feynman diagram for VEPP life cycle. (Time advances in the upward di- 
rection in this figure^.) The bottom vertex is spacetime point xj, the top vertex is 
spacetime point X 2 . The maximum extent of this process along the time axis is t// e . 

(b) Semiclassical VEPP life cycle model. The magnetic moments fj, e ± are assumed to 
point along the axis of separation. The time ti corresponds to \t ve in the case of a 
VEPP. 

Figure 2. Magnetic field lines of a positron at the origin. The field is symmetric with 
respect to rotation about the vertical axis. An example location of a virtual electron is 
indicated by the vertical arrow, representing its magnetic moment n e - , which must al- 
ways be parallel with the magnetic moment of the positron because angular momentum 
is conserved. The dashed lines, which become a pair of cones when rotated about the 
vertical axis, define the volume where energy conservation permits virtual electrons to 
exist: within the cones electrons may move, but outside the cones this violates the en- 
ergy equation (2). 
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TABLE I. 

Elementary Particles and Their VAP Properties 


Particle 

Mass (g) 

qje 


■^moi (cm) 

tv (s) a 

tH (s) 




A 




Electron 

9.11 x 10“ 28 

-1 


4.7 x 10~ 12 

2.7 x lO" 22 

3.2 x lO' 22 

Muon 

1.88 x 10" 25 

-1 


2.4 x 10" 14 

1.3 x 10~ 24 

1.55 x 10“ 24 

Tau lepton 

3.16 x 10" 24 

-1 


1.3 x 10~ 15 

7.8 x 10“ 26 

9.23 x 10" 26 

Down quark 

5.9 x 10- 25 

1 

3 

0.52 

2.1 x 10' 14 © 

1.4 x 10~ 24 © 





0.0048 

2.1 x 10" 14 

1.4 x 10~ 24 

5.0 x lO" 25 

Up quark 

5.9 x IQ" 25 

2 

3 

0.52 

2.1 x 10~ 14 © 

© 

CN 

1 

o 

X 





0.0193 

2.1 x 10" 14 

1.4 x 10“ 24 

5.0 x 10" 25 

Strange quark 8.9 x 10 25 

1 

3 

0.38 

1.0 x 10~ 14 © 

6.7 x 10" 25 © 





0.0123 

1.0 x 10" 14 

6.7 x 10" 25 

3.3 x lO’ 25 

Charm quark 

2.67 x 10“ 24 

2 

3 

0.22 

1.9 x 10“ 15 © 

1.3 x 10~ 25 © 





0.26 

2.3 x 10“ 15 

1.5 x lO' 25 

1.1 x lO" 25 

Bottom quark 8.4 x 10 24 

1 

3 

0.19 

5.3 x 10" 16 © 

3.6 x 10“ 26 © 





0.098 

5.8 x 10" 16 

3.7 x 10" 26 

3.5 x 10" 26 

Top quark 

3.1 x 10" 22 

2 

3 

0.11 

8.4 x 10“ 18 © 

5.6 x 10~ 28 © 





2.0 

1.4 x lO" 17 

8.6 x lO” 28 

9.4 x lO" 28 

W boson 

1.42 x 10- 22 

l 


3.0 x 10“ 17 

1.7 x 10- 27 

2.1 x 10" 27 

Z boson 

1.62 x 10" 22 

0 


2.6 x lO -17 

1.5 x 10" 27 

1.8 x 10- 27 


° Two values are given for quarks, the upper number marked with the © symbol to indi- 


cate that the color- Coulomb potential only was used to derive the value; the lower num- 
ber is found by including both color-Coulomb and spin-spin potentials. 



Figure 1(a), The Lifetimes of Virtual 
Antiparticle Pairs L 
by David Batchelor 
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Figure 1(b), The Lifetimes of Virtual 
Antiparticle Pairs I 
by David Batchelor 
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Figure 2, Semiclassical Models for Virtual 
Antiparticle Pairs, the Unit of Charge e and 
the QCD Coupling a,, by David Batchelor 
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